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Wetting failure is of vital importance to many physical phenomena, such as industrial coating and drop 
emission. Here we show when and how the surface roughness promotes the destabilization of a moving 
contact line on a hydrophilic surface. Beyond the balance of the driving force and viscous resistance where a 
stable wetting interface is sustained, wetting failure occurs and is modified by the roughness of the surface. 
The promoting effect arises only when the wetting velocity is high enough to create a gas-liquid-solid 
composite interface in the vicinity of the moving contact line, and it is a function of the intrinsic contact 
angle and proportion of solid tops. We propose a model to explain splashes of rough solid spheres impacting 
into liquids. It reveals a novel concept that dynamic wetting on hydrophilic rough surfaces can be similar to 
that on hydrophobic surfaces, and brings a new way to design surfaces with specific wetting properties. 

Wetting phenomena are ubiquitous in biological activities 1 , liquid coating 2 " 6 , microfluidics 7 , etc. As a 
limiting case, wetting failure gives rise to abundant physical problems such as solid-liquid impact 
splashes 8 ' 9 , drop emission 10 , and air entrainment in industrial coating 2 " 6 . Wetting dynamics on rough 
surfaces dates back to the experimental work by Cazabat 11,12 . He observed the spreading of the silicone oil on both 
smooth and rough glass surfaces, and revealed complex spreading regimes in addition to the traditional regimes of 
spreading driven by capillarity and gravity. With further experimental and theoretical studies of the contact line 
dynamics on ideally smooth surfaces 13 " 19 , the influences of roughness on dynamic wetting have been explored 
recently, e.g., hemi-wicking of superhydrophilic spreading 20 ' 21 , coating limits on rough surfaces 5 ' 6 , pinning poten- 
tial by defects 22 , and topographically driven spreading 23 " 25 . However, most present studies only consider one 
certain aspect such as superwetting limit 20,22 " 25 , superhydrophobic limit 26 , and low capillary number limit 20,22 " 25 . 
Clarke 5 and Benkreira 6 carried on experiments to determine how the roughness modifies the wetting failure 
systematically outside of those limits; however, no explanation was provided to describe when and how the 
roughness effects occur on an arbitrary hydrophilic surface. 

The impact of a solid sphere onto a liquid pool provides a very ingenious and economical way for the study 8 on 
the effect of surface roughness on dynamic wetting. The promotion of the wetting failure is observed in our 
experiments and there are currently no convincing explanations available. In fact, the enhancement of dynamic 
wetting has been observed in previous studies 23 " 25 , which showed dynamic influences of roughness on the 
superhydrophilic drop spreading and discovered a simple modified power law relationship between the edge 
speed and dynamic contact angle. The experimental and theoretical work, albeit falling into the limits of low 
spreading speed and superwetting regimes, indicates that the wetting dynamics is strongly dependent on the 
surface roughness. However, when the static contact angle is finite or the wetting velocity is high (corresponding 
to a relatively high capillary number), what will happen? Will the wetting process be enhanced or inhibited by the 
surface roughness? 

In this paper, by observing the spreading film generated by the solid-water impact experimentally, similar to 
what has been performed by other groups before 8,9,27 , we are going to explore when the promoting effect on 
wetting failure will occur by providing a criterion to predict the appearance of a liquid-gas-solid composite 
interface in the vicinity of the moving contact line (MCL). We will also show how the promoting effect works 
based on the contact-line instability induced by the imbalance between the capillary driving force and viscous 
friction force 8,10,17,19,24 . 



Results 

The stability of the thin liquid film climbing on an impact body, which is relevant to surface roughness, deter- 
mines whether an air cavity will be created or not 5,6,19 (as illustrated in Fig. 1). On impacting into a water pool, a 
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Figure 1 | Splashes generated by a partly roughened glass sphere impacting into water pool with velocity 3.6 ms" 1 . When a glass sphere partly 
roughened with sandpapers impacts into the initially quiescent water pool (a), a water film generated by the impact propagates in different manners. On 
the roughened patch, the film detaches from the sphere surface and induces a big cavity, as illustrated by the right half of the sphere in (b); however, on the 
smooth periphery, the film moves well up conforming to the sphere surface, as illustrated by the left half of the sphere in (b) . The behavior of the water film 
is a dynamic wetting phenomenon and the static wetting property of the surface is a key factor in determining the stability of MCL 81018 marked in (a). 
Static contact angles of the smooth and rough part of the sphere are about 49.5° and 38.2°, respectively; thus, the roughness promotes wetting failure 
significantly in comparison with previous predictions and experiments 81019 , (b) is the schematic of splashes on a half roughened surface. V 0 and Vf are the 
impact velocity of sphere and wetting velocity of the water film, respectively. V}-~ £ V 0 8,34 ' 35 . The thickness of splash, 5 film , is on the scale of a millimeter 35 . 



smooth sphere partly roughened by the abrasive sanding process 
shows complex splash patterns: a distinct air cavity is created along 
the sanding patch where a typical sawtooth-shaped MCL is always 
observed, while a well- developed MCL in the smooth portion moves 
smoothly up along the sphere periphery. Since the physical and 
chemical parameters in and out of the patch are identical, except 
for the static contact angle and surface roughness, it can be sure that 
these two parameters make the difference. Detailed results on the 
wetting failure of statically hydrophilic spheres with various rough- 
ness and materials are depicted in Fig. 2 (see Tables SI -4 for more 
details), which indicates a strong dependence of wetting failure on 
the surface roughness. 

Discussion 

Figure 2a demonstrates that the surface roughness plays more evid- 
ent roles than the static contact angle. A crucial question is evoked: 
how could the surface roughness decrease the wetting failure velocity 
in such a generic manner? Multi- scale forces are involved in prop- 
agating the water film 8,1019,24 : external force, capillarity, and viscous 
stresses. In the inner region as depicted in Fig. 2c, the viscous friction 
force and capillarity dominate the bending of the interface. While in 
the outer region out of a characteristic length h c from the surface wall, 
capillarity and the external forces, such as inertia and gravity, are 
dominant ones to sustain the balance. The characteristic length, h a 
dividing the inner and outer region, is on the scale of the capillary 

length \/y/pg, where y is the surface tension between the liquid and 
air, p is the density, and g is the gravitational acceleration. A global 
stability of the advancing water film is achieved by an intimate 
matching of the boundary restrictions and outer interface through 
the inner region. The outer interface can always bend its profile to 
hold the outer balance, but there is a limitation for the inner region to 
achieve the matching 17 . Therefore the wetting failure always starts 
from the inner one. Consequently, to characterize the dynamic wet- 
ting failure, we cast our attention on the inner region. 

In the inner region 81019 ' 24 , the capillary driving force could be 
considered as the rate of change of surface free energy. Based on that, 
we will show how the dynamic wetting behavior is modified by the 
topography into two regimes. 

• When a liquid film creeps on a rough surface at a low velocity, it 
advances conforming to the surface topography 5 . Thus, the 



change of surface free energy, when the contact line propagates 
by a distance dz, is modified as dE = r(y S L ~ ysv)dz + ycosO d dz. In 
this equation, r, y SL , and y sv are the ratio of the actual surface to 
the projective area, the surface tension between the solid and the 
liquid, and the surface tension between the solid and the air, 
respectively. By employing the Wenzel equation 21 , we deduce 
the driving force dE/dz = y(cosO w — cosO d ). On a smooth surface 
with a static contact angle 6 Y , this driving force is simply y(cosO Y 
— cosO d ). As 0 W is smaller than 6 Y , the driving force is enhanced 
by the surface roughness. More importantly, if the variation of the 
viscous dissipation due to the surface roughness is ignored, the 
dynamic wetting behavior is driven by the topography 24 . 
• However, when the liquid film skips from peak to peak on rough 
surfaces, composite liquid- air- solid interfaces will be induced by 
the MCL 5 . As the MCL moves by dz, the surface free energy 
changes by dE =f s (ysL ~ ysv)dz + (1 —f s )y dz + ycosO d dz, where 
f s is the fraction of solid tops. Consequently, the capillary driving 
force in this case is y(cos9 C B ~ cosO^), where 0 C b is calculated by 
the Cassie-Baxter equation 21 . If the viscous dissipation change is 
ignored, the dynamic wetting process is inhibited because the 
driving force is decreased. 

Since the driving force in the inner region can always be written as 
y (cosO m — cosO d ), here we provide a tentative explanation for the 
ubiquitous microscopic contact angle 0 m in regard to the surface 
roughness. 0 m depends on the instantaneous state of the MCL: if 
the MCL propagates clinging to the topography, 0 m can be obtained 
by Wenzel equation; if the MCL propagates on tops of the topo- 
graphy, 0 m is calculated by Cassie-Baxter equation. The understand- 
ing is crucial to filling the full scope of wetting dynamics, and brings 
out a subtle mechanism in which the surface heterogeneity may 
hinder or facilitate the wetting motion. Before going further, another 
two significant questions arise: when will the effect of promoting and 
postponing the wetting failure occur, and what is the effect of surface 
roughness on the viscous friction force? 

To obtain a reasonable driving force on rough surfaces, we approx- 
imate the rough surface as a smooth surface sculptured with arrays of 
vertical capillaries (as shown in Fig. 2b and c) and analyze the cri- 
terion at which the Wenzel to Cassie-Baxter state transition appears. 
Therefore, it is worth considering the case that a liquid film passes 
over one single capillary and comparing two time scales: t c which is 
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Figure 2 | Wetting failure velocity versus static contact angle and surface roughness (a) and schematic of a water film moving on the surface 
topography (b), (c). The blue curve in (a) is extracted from the experimental data and theory 8 ' 10 . Inset images are experimental snapshots of 30 ms after 
impact for two identical aluminum spheres, with diameter 25.4 mm and impact velocity 2.9 ms -1 , except for their surface roughness (Left, R z = 2.37 (am; 
Right, R z = 9.35 (am). Above the blue curve in (a) 810,19 , a big cavity is always observed for a smooth surface (the cavity forming entry), while below the 
curve, no cavity is created (the silent entry). In sharp contrast to this prediction, the surface roughness on a statically hydrophilic surface decreases the 
wetting failure velocity, V* y significantly, as depicted in the inset of (a). Further, a rough surface is artificially approximated as a smooth surface patterned 
by capillaries (b), and the characteristic depth h and diameter L of capillaries are determined experimentally and statistically. Near the moving contact line 
(b), two regions 10 ' 19 are defined (c) according to the transition length h c . In the outer region, the profile of the film is dominated by external forces and 
capillarity, and the dynamic contact angle 6 d is observed 1016 " 19 . In the inner region, the interface is mainly portrayed by the capillarity and viscous friction 
force, and it meets the solid at a truncated length scale h m with a microscopic contact angle 0 m 10 ' 16 " 19 . When the water film moves on the topography, not 
only will the liquid front advance in the direction parallel to the surface with velocity Vp but also it will fill up the valleys of the topography in the vertical 
direction with a speed V c . 



the time scale of filling up the gap with the liquid and Xf which is the 
time scale of the MCL moving over the opening of the capillary with 
the advancing velocity of Vf. The ratio, E = xjxj > 1, is therefore 
proposed as the criterion for establishing the forced Cassie-Baxter 
state, in which case air is sealed in topography valleys and thus the 
wetting film front moves on a liquid-air-solid composite surface. The 
invading of the liquid into the valley is considered as a spontaneous 
wetting process in a capillary, which is shown in Fig. 3. Since x c is 
calculated based on a free capillary rise 28 , it is rather under- estimated 
due to the neglect of the impedance of the contact area, inlet sharp 
edge, air pressure, etc. As to the viscous friction force, the influence of 
topography is expected to be marginal 8 ' 29 ' 31 " 33 . By taking the two 
aspects, i.e. capillary driving force and viscous resistance, into con- 
sideration, the dynamic contact angle on the rough surface is larger 
than that on the smooth one at the same wetting velocity if the two 
surfaces (smooth and rough) have the same static contact angle 0 e , 
when the Wenzel to Cassie-Baxter transition occurs. As a result, the 
dynamic wetting is inhibited. 

Now let us go back to the results summarized in Fig. 2a. First of all, 
for spheres with almost ideally smooth surfaces (the smooth glass 
sphere), the critical velocity agrees well with the theoretical predic- 
tion 810 , where the driving force is derived based on the deviation of 
the dynamic contact angle from static one and the viscous resistance 
including the contributions of the liquid (F L ~ Cfi L Vf) and the air (F air 
~ (3fi air l/[n — O d ])vf) 8,29,30 . However, as to rough hydrophilic spheres, 
the critical velocity decreases drastically with the increasing of rough- 
ness plotted in Fig. 2(a), implying that effective driving force is 
reduced or the overall viscous friction force is increased in the dia- 
gram of hydrodynamic model. To find out the reasons, we detected 
sphere surfaces with roughmeter and SEM, and reconstructed their 
3-D topography with the Laser Scanning Confocal Microscope. 



Results show that E > 1 (See Supplementary Information and 
Table S5 for more details) when the wetting velocity approaches 
the critical velocity. This implies that the MCL moves on an air- 
liquid-solid composite surface depicted in Fig. 2b. In this situation, 
the capillary driving force is decreased from y (cos0 e — cos0 d ) to 
y(cos0 C B ~ cos0 d ), as aforementioned. On the other hand, following 
Ybert's procedure 32 , the slip length induced by the surface roughness 

z k L Regime I: pure acceleration regime, z~t 2 
Regime II: linear regime, z~t 




Figure 3 | Capillary tube rise of low viscosity liquid 27 , r, z y p, y, rj y and V c 

are the radius of a capillary tube, height of meniscus, liquid density, surface 
tension, viscosity, and meniscus rising velocity, respectively. The liquid-air 
interface of the liquid pool is located at z = 0 and the liquid starts to intrude 
into the capillary at t = 0. If the liquid is more viscous than gp 3/2 r 5/2 y~ 1/2 , the 
capillary tube imbibition process can be divided into three stages 28 : 
acceleration, linear, and Washburn regimes (see Supplementary 
Information for the description of the capillary rise). Accordingly, E for a 
given system can be calculated (Table S5). 
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Figure 4 | Air entrainment velocity after revising the driving force vs. 
theoretical predition. The x-axis coresponds to the Cassie-Baxter contact 
angle; the y-axis coresponds to the range of air entraiment velocities. Blue 
line is the theoretical prediction 810 . Different symbols stand for differet 
surfaces: smooth glass sphere (diameter = 48 mm, R z = 0.96 (am; cyan 
circle), rough glass sphere( diameter = 48 mm, R z = 8.83 (am; cyan 
square), aluminum sphere I (diameter = 25.4 mm, R z = 2.37 (am; black 
circle), aluminum sphere II (diameter = 25.4 mm, R z = 2.83 (am; black 
diamond), aluminum sphere III (diameter = 25.4 mm, R z = 5.61 (am; 
black triangle), aluminum sphere IV (diameter = 25.4 mm, R z = 9.16 (am; 
black square), aluminum sphere V (diameter = 30 mm, R z = 1.76 (am; 
green circle), aluminum sphere VI(diameter = 30 mm, R z = 2.3 (am; 
green diamond), aluminum sphere VII(diameter = 30 mm, R z = 
3.25 (am; green triangle), aluminum sphere VIII (diameter = 30 mm, R z = 
7.02 urn; green square), rough nylon sphere (diameter = 50 mm, R z = 
21.77 (am; red square). Intrintic contact angle 0 y are experimentally 
extracted as 37.4°, 76.5°, 59.2°and 71.2° for glass sphere, aluminum sphere 
I-IV (diameter = 25.4 mm), aluminum sphere V-VII (diameter = 
30 mm) and nylon sphere, respectively. 

is no more than 5 \im, indicating that overall viscous resistance will 
not be altered more than 5% 29 " 33 . Thus the viscous dissipation is 
treated similar to that on a smooth surface or a hydrophobic surface 
with micro-roughness 8 ' 29 . An obvious drop of the critical velocity is 
then expected. For rough surfaces, Duez and Ybert's threshold of 
impact velocity 81019 is explicitly modified as 
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where C, l,g 0 are numerical pre-factors introduced to account for the 
liquid, air, and overall viscous dissipations 8,29,30 . These factors can be 
roughly determined through other more fundamental experiments. 
The value of / falls in the range of about 15 ~ 20 8,29,30 , depending on 
details of interactions of the liquid and the solid substrate. g 0 and C 
can be determined by balancing the viscous friction force and the 
capillary driving force when / is fixed. The dependence ofg 0 and C on 
/ is numerically presented in supplementary information; their 
experimental values fall in the expected range 815 ' 29 ' 30 , verifying the 
validity of the model, q is the ratio of the splash film velocity and 
impact velocity (~2 8,9,34,35 ). £ is an experimentally and theoretically 
determined constant in that the increase of the capillary driving force 
with the surface hydrophilicity is almost neutralized by a simulta- 
neously enhanced viscous resistance 8,19,29 for hydrophilic surfaces. 
From equation (1), the critical capillary number remains the same 



for the revised contact angle 0 C b for an arbitrary rough hydrophilic 
surface. It agrees with the previous theoretical and experimental 
conclusions 8,10,19 . As is shown in Fig. 4 (Calculation parameters can 
be found in Tables S5 and Supplementary Information), the modi- 
fied theoretical prediction agrees well with experimental results, 
which validates the feasibility of our model. 

In conclusion, we have demonstrated a hidden mechanism for the 
wetting failure promoted by surface roughness. The surface rough- 
ness can modify wetting dynamics by changing the capillary driving 
force in the inner region. At a low wetting velocity, the roughness of a 
hydrophilic surface enhances the wetting process. As the wetting 
velocity increases to a critical value, the Wenzel to Cassie-Baxter 
transition occurs and the surface roughness starts to inhibit wetting. 
This mechanism is different from the conventional theory on rough- 
ness-induced wettability, and it provides a new prospective to under- 
standing the ubiquitous micro -contact angle and opens a novel way 
to design specific rough surfaces. 



Methods 

Experimental setup. The apparatus used for determining the critical velocity of the 
water impact is detailed in our previous work 27 . High speed camera was set at the 
frequency of 2000 frames per second, and the shot images were used to identify 
wetting failure transition and calculate the impacting velocity. The critical velocity for 
each sphere was obtained by averaging high speed camera data from at least 3 runs to 
ensure the repeatability of results. Before each run, spheres were cleaned in the 
deionised water and dried naturally in lab conditions. Experimental temperature was 
stabilized at 20 ± 2°C. 

Surface treatment and measurement. Spheres used were aluminum (diameter, 
25.4 mm and 30 mm), nylon (diameter, 50 mm) and glass spheres (diameter, 
48 mm). Different addictive materials were added to aluminum spheres with 
different diameters. The surface roughness was acquired by an abrasive sanding 
process with the grit size in the range of 150 ~ 2000-grit and quantitatively measured 
by a roughmeter. Measurement parameters of the roughmeter were set as L t = 
1.75 mm, L c = 0.25 mm, and n = 5. Static contact angles of different surfaces were 
measured by a contact angle meter, OCA15EC of DataPhysics Instruments, using the 
standard sessile drop method. Surfaces were detected further by SEM and Laser 
Scanning Confocal Microscope, Olympus LEXT OLS4000, and their 3-D patterns 
were reconstructed according to the data obtained. The material ratio rate with regard 
to the altitude of the 3-D surface was used to determine the ideal smooth surface (the 
reference plane), and then the surface characteristic length and the ratio of solid tops 
were calculated based on the reference plane. The length scale of the detection area 
was at least five times of the surface characteristic length. 
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